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Abstract 

Pseudoscalar meson decay leads to an entangled state of charged leptons (fi, e) and massive 
neutrinos. Tracing out the neutrino degrees of freedom leads to a reduced density matrix for the 
charged leptons, whose off-diagonal elements reveal charged lepton oscillations. Although these 
oscillations decohere on unobservably small time scales < 10~ 23 s, they indicate charged lepton 
mixing as a result of common intermediate states. The charged lepton self energy up to one 
loop features flavor off-diagonal terms responsible for charged lepton mixing: a dominant "short 
distance" contribution with W bosons and massive neutrinos in the intermediate state, and a 
subdominant "large distance" contribution with pseudoscalar mesons and massive neutrinos in the 
intermediate state. The mixed u — e propagator cannot be completely diagonalized by a unitary (or 
bi unitary) transformation as a consequence of the different spinor structure between the kinetic 
and mass terms. Mixing angle(s) are GIM suppressed and are momentum and chirality dependent. 
The negative chirality mixing angle near the muon mass shell is 9l(M^) oc Gf^2 ^^j^^je w ith 
mj the mass of the neutrino in the intermediate state. Recent results from TRIUMF suggest 

an upper bound 9l(M^) < 10 -14 ( Ms/100 MeV\ for one generation of a heavy sterile neutrino 



with mass M$- We obtain the wavefunctions for the propagating modes and discuss the relation 
between the lepton flavor violating process \x — > ej with charged lepton mixing, highlighting that 
a measurement of such process implies a mixed propagator u, e and suggest further contributions 
to this process as a consequence of mixing with momentum dependent mixing angles. 
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I. INTRODUCTION 



Neutrino masses, mixing and oscillations are the clearest evidence yet of physics beyond 
the standard model Oscillations among three "active" neutrinos with 5m 2 = 10~ 4 — 

10~ 3 eV 2 , for atmospheric and solar oscillations respectively, have been firmly confirmed 
experimentally (see the reviews [H[] - fill] ) - 

However, several experimental hints have been accumulating that cannot be interpreted 
as mixing and oscillations among three "active" neutrinos with Sm 2 ~ 10~ 4 — 1CT 3 . Early 
results from the LSND experiment (l2| have recently been confirmed by MiniBooNE run- 
ning in antineutrino mode[13|], both suggesting the possibility of new "sterile" neutrinos 
with 5m 2 ~ eV 2 . The latest report from the MiniBooNE collaboration on the com- 
bined Vp — > v e and — > V e appearance data is consistent with neutrino oscillations with 
0.01 < Am 2 < 1.0 eV 2 . This is consistent with the evidence from LSND antineutrino 
oscillations (l2l|. which bolsters the case for the existence of sterile neutrinos; however, com- 
bined MiniBooNE/SciBooNE analysis [15] of the I7 M disappearance data are consistent with 
no short baseline disappearance of V^. Recently, a re-examination of the antineutrino flux 16[ 
in anticipation of the Double Chooz reactor experiment resulted in a small increase in the 
flux of about 3.5% for reactor experiments leading to a larger deficit of 5.7% suggesting 
a reactor anomaly [l7| ■ If this deficit is the result of neutrino mixing and oscillation with 
baselines L < 10 — 100 m, it requires the existence of at least one sterile neutrino with 
8m 2 > 1.5 eV 2 and mixing amplitude sin 2 (2#) ~ 0.11513]. Taken together these results 
may be explained by models that incorporate one or more sterile neutrinos that mix with 

although 



the active ones[18l]-[25j including perhaps non-standard interactions 26]; although, there is 
some tension in the sterile neutrino interpretation of short-baseline anomalies [2 71] . A com- 
prehensive review of short baseline oscillation experiments summarizes their interpretation 
in terms of one or more generations of sterile neutrinos 28 . 



Recently it has been pointed out that the presence of sterile neutrinos may induce a 
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modification of the recently measured angle #13 [21 

Furthermore, the latest analysis of the cosmic microwave background anisotropies by 
WMAP[31] suggests that the effective number of neutrino species is N e ff = 4.34 ± 0.86 and 
^2(m u ) < 0.58 eV^, bolstering the case for sterile neutrino(s) with m < eV. Complementary 
cosmological data suggests that N e ff > 3 at the 95% confidence level 3^]; although, accom- 
modating an eV sterile neutrino requires a reassessment of other cosmological parameters [33|. 
For recent reviews on "light" sterile neutrinos see ref. 34J. Sterile neutrinos with masses in 
the ~ keV range may also be suitable warm dark matter candidates [351]- [40j and appealing 
models of sterile neutrinos provide tantalizing mechanisms for baryogenesis[4l|. 

These hints motivate several experimental proposals to search for sterile neutrinos (see 
the reviews in ref. [34|). Various experimental searches have been proposed, such as Higgs 
decay and matter interactions of relic sterile neutrinos the end point of /3-decay in 
187 Re with a value of Q = 2.5keV[H, 0], electron capture decays of 163 Ho -)- 163 .Ch/jll] 
and 8 Li production and decayj46|. More recently, the focus has turned on the possible 
new facilities at the "intensity frontier" such as project X at Fermilab j47| . alternative high 
intensity sources [H, 48] and recent proposals to study sterile-active oscillations with pion 



and kaon decay at rest (DAR)[49j, [50| or muons from a storage ring 51] as well as the 



possibility of discrimination between heavy Dirac and Majorana sterile neutrinos via \AL 
2 processes in high luminosity experiments [521] . which is summarized in a recent review 28 
Although the recently reported analysis of the phase II data of the Mainz Neutrino Mass 
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Experiment [53] found no evidence for a fourth neutrino state tightening the limits on the 
mass and mixing of a fourth sterile species, the possibility of a heavy sterile species is 



still actively explored [54], [55j. More recently the PIENU collaboration at TRIUMF[56] has 
reported an upper limit on the neutrino mixing matrix element \U e i\ 2 < 10~ 8 (90% C.L.) in 
the neutrino mass region 60 — 129MeV/c 2 . 

In this article we focus on complementary consequences of sterile neutrinos in the form 
of charged lepton mixing phenomena. The discussion of whether or not charged leptons 
oscillate has been controversi al [oTI] - (63| and, more recently, this question was addressed from 
the point of view of coherence [641] which highlights that, while oscillations are possible, they 
lead to rapid decoherence and no observable effects. Muon-antimuon oscillations via massive 
Majorana neutrinos have been studied in ref.(65[; however, to the best or our knowledge, 
the issue of charged lepton (fi — e) mixing (we emphasize mixing over oscillations) has not 
yet received the same level of attention. Although in ref. 66|| charged lepton mixing and 



oscillations as a consequence of neutrino mixing was studied in early Universe cosmology at 
temperatures C T < Myy where it was argued that medium effects enhance charged 
lepton mixing, the question of charged lepton mixing in vacuum and as a consequence of 
possible new generations of sterile neutrinos has not yet been discussed in the literature and 
is the main motivation of this article. 

Furthermore we discuss the relationship between the lepton flavor violating decay /i — > 
and charged lepton mixing in terms of self-energies and propagators that mix // and e. 
Charged lepton violation is the focus of current experimental searches [671. 68] and a recent 



experimental proposal [69] to search for charged lepton flavor violation via the coherent 
conversion process \x — N — > e — iV at Fermilab. 

Goals: In this article we study both charged lepton oscillations and mixing as a con- 
sequence of intermediate states of mixed massive neutrinos and discuss the relationship 
between charged lepton mixing and charged lepton flavor violating processes. 

• a) Oscillations: In a recent article[70j (see also 71, 72||), we have provided a non- 
perturbative quantum field theoretical generalization of the Weisskopf-Wigner method 
to understand the correlated quantum state of charged leptons and neutrinos that con- 
sistently describes pion/kaon decay in real time. Knowledge of this state allows us to 
obtain the reduced density matrix for charged leptons by tracing out the neutrino de- 
grees of freedom. The off diagonal density matrix elements in the flavor basis contains 
all the information on charged lepton (/i, e) coherence and oscillations. 

• b) Mixing: Charged lepton oscillations evidenced in the reduced density matrix are a 
consequence of a common set of intermediate states that couple to the charged leptons. 
We then study the charged current contribution to the one-loop self-energy which 
couples charged leptons to an intermediate state of mixed massive neutrinos. The 
self-energy unambiguously determines the propagating states and explicitly describe 
charged lepton mixing. We obtain the mixed propagator while extracting the mixing 
angles and analyze the propagating modes and their wavef unctions. These results 
motivate us to address the relation between lepton flavor violating transitions such as 
(/, —> ery and charged lepton mixing. 



Brief summary of results: 
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a:) The quantum state of charged leptons and neutrinos from (light) pseudoscalar 
decay is a correlated entangled state from which we construct the corresponding (pure 
state) density matrix. Under the condition that neutrinos are not observed, we trace 
over their degrees of freedom leading to a reduced density matrix for the charged 
leptons. Because we focus solely on decay of ir,K the reduced density matrix de- 
scribes the remaining /j,, e. Integrating out the unobserved neutrinos leads to a re- 
duced density matrix that is off-diagonal in the flavor basis. The off-diagonal matrix 
elements describe charged lepton mixing and exhibit oscillations with typical frequency 
E/j.(k) — E e (k) > 0{m ll — m e ) ~ m M ~ 1.6 x lO 23 ^^ 1 which are unobservable over any 
experimentally relevant time scale and lead to rapid decoherence. This conclusion 



agrees with a similar observation in ref . 64] . While these fast oscillations lead to deco 



herence over microscopic time scales, we recognize that the origin of these oscillations 
are a common set of intermediate states akin to neutral meson oscillations. 

Recognizing that the origin of oscillations are intermediate states that are common 
to both charged leptons, we obtain the self-energy contributions and the full mixed 
propagator for the /i, e system. Mixing is a direct result of charged current interac- 
tions with intermediate neutrino mass eigenstates. As in the case of neutral meson 
mixing, we identify "short" and "long" distance contributions to the flavor off-diagonal 
self-energies. The dominant "short" distance contribution corresponds to the interme- 
diate state of a W ± and neutrino mass eigenstates, whereas the subdominant "long" 
distance contribution is described by an intermediate state of ir, K and a neutrino 
mass eigenstate. We calculate explicitly the short distance and provide an estimate 
for the long distance contributions. Unitarity of the neutrino mixing matrix entails 
a Glashow-Ilioupoulos-Maiani (GIM) type mechanism that suppresses charged lepton 
mixing for light or nearly degenerate neutrinos, thus favoring heavy sterile neutrinos 
as intermediate states. 

We obtain the flavor off-diagonal charged lepton propagator and analyze in detail 
the propagating modes, fi — e mixing cannot be described solely in terms of a 
local off-diagonal mass matrix but also off-diagonal kinetic terms which are four- 
momentum dependent and contribute to off-shell processes. Mixing angles are GIM 
suppressed and both chirality and four momentum dependent. The largest angle cor- 
responds to the negative chirality component near the muon mass shell, it is given 
by 9l(M^) oc GF^jUpjUj e rrij where rrij is the mass of the intermediate neutrino. 
Therefore charged lepton mixing is dominated by intermediate states with mixed 
heavy neutrinos. Assuming one generation of a he avy sterile neutrino with mass 
Ms and extrapolating recent results from TRIUMF[56] we obtain an upper bound 

9 L (M^) < 1CT 14 ^Ms/100 Me . We obtain the propagating eigenstates of charged 

leptons and show that the full propagator cannot be diagonalized by a simple unitary 
(or bi-unitary) transformation. 

The relationship between charged lepton mixing and the lepton flavor violating decay 
fi — > e7 is discussed in terms of the mixed charged lepton self-energies and possible 
observational effects in the form of further contributions to [l — > cy are discussed. 
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II. REDUCED DENSITY MATRIX: CHARGED LEPTON OSCILLATIONS 

In ref. [70[ the quantum field theoretical Weisskopf-Wigner (non-perturbative) method has 
been implemented to obtain the quantum state resulting from the decay of a pseudoscalar 
meson M, (pion or kaon). It is found that such state is given by (see [70] for details and 
conventions), 

|Mr(t))) = e-^'^e-^'^lMziO)) - (II.l) 
where \^i^(t)) is the entangled state of charged leptons and neutrinos given by 



*»,*(*))= \u a j^ajM aj (k,q,h,h') ^[k, q;t\ 

q,aj,h,h' v 



x e- i( ^ (fc)+ ^ (9)) *|C(^fc))|^(^,-g)) \ ; k = p + q, (11.2) 

where 



JF aj [q,p,h, h';t] 



1 _ e -i(^ f (p)-S Q (fc)-^(g)-i^-)t 
E^ip) - E a (k) - Ej{q) - i^- 



(II.3) 



where A4 a j(k, q, h, h'), IT Qj j(g, k) are the production matrix elements and phase space factors 
respectively, 

Maj(k, q, h, ti) = F M Uc h h$)'fW j , h .(q)p l t (II.4) 

U a Aq, k) = 1 (11.5) 

^J&V E M (p)E a (k)Ej(q) 

In these expressions F M is the pion or kaon decay constant, U.;V^h'(q) are the spinors 
corresponding to the charged lepton a and the neutrino mass eigenstate j (for notation and 
details see ref. [7(|). The leptonic density matrix that describes the pure quantum entangled 
state of neutrinos and charged leptons is given by 

P iAt) = \^iAt))(^iAt)\ (n.6) 

If the neutrinos are not observed, then their degrees of freedom must be traced out in the 
density matrix and the resulting density matrix is no longer a pure state, 

Pf(t) = Tr vPliU (t) (H.7) 

Considering only light pseudoscalar decay tt, K, the only charged leptons available are 
H, e. For a fixed helicity h and momentum k of the charged leptons, the reduced density 
matrix is given by 

pf(t) = p ee (h,k,t)\e^) (e- z \+Pw(h,k,t)\fi- n ) (ji~ n \ (H.8) 



+ p e ,{h,k,t)\e-J {^\e*^-^ + p^(hXt)\^ 



> e -J e -i(fi M (fe)--B e (fe))t 
v h,k 
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The diagonal density matrix elements in the //, e basis describe the population of the 
produced charged leptons whereas the off-diagonal elements describe the coherences. The 
diagonal matrix elements p aa , a = /i, e are given by 

pL(t) = J2\ u ^\ 2BR ^^A 1 - e ' rMt ] ; « = ^ e ( n - 9 ) 

where BR are the branching ratios T M ^ ia U ./T M and we have used some results obtained 



in ref.[70|. The off diagonal elements do not have a simple expression, however the most 



important aspect for the discussion is that these density matrix elements are of the form 

p% « E M ^ M li ; 4 = W > ( IL1 °) 

5 

namely, they describe the process M — > auj followed by a "recombination" -type process 
Mvj — > (3 thereby suggesting the intermediate state Muj common to both matrix elements. 
For r M i > 1 the reduced density matrix (for fixed h, k) in the charged lepton basis is of the 
form 



P R 



A ee A ee -K^(k)-E a (k))t 

A e i{E^(k)-E e {k))t a 
^e/i MA 1 



(11.11) 



This tells us that there will be fi <^ e oscillations. However, these oscillations occur 
with large frequencies E^{k) — E e (k) > 0{m^ — m e ) ~ m M ~ 1.6 x 10 23 s~ 1 and are unob- 
servable over any experimentally relevant time scale. This conclusion agrees with a similar 
observation in ref . [641] . 

Although these oscillations average out over relevant time scales and are experimentally 
unobservable, an important issue is their origin. The mixing between charged leptons arises 
from the fact that they share common intermediate states. In the case studied above, 
the common intermediate state corresponds to a pseudoscalar meson and a neutrino mass 
eigenstate. 

From this point of view, the physical origin of the oscillations is found in mixing of the 
charged leptons from the fact that they share common intermediate states. This is in fact 
similar to the oscillations and mixing in the KqKq system. The obvious difference with 
this system is that, in absence of weak interactions, Kq and Kq are degenerate and this 
degeneracy is lifted by the coupling to the (common) intermediate states, thereby leading 
to oscillations on long time scales. 

The conclusion of this discussion is that charged lepton oscillations are a result of their 
mixing via a set of common intermediate states. The off-diagonal density matrix elements are 
the lowest order in a perturbative expansion, therefore they do not reveal the full structure 
of the mixing phenomenon. 

If charged leptons mix via a common set of intermediate states, the correct propagating 
degrees of freedom correspond to the full charged lepton propagator which requires the self- 
energy correction. Such self-energy will reflect the mixing through the intermediate states. 
Whereas oscillations average out the off-diagonal density matrix on short time scales, the 
main physical phenomenon of mixing is manifest in the description of the true propagating 
modes. 
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III. CHARGED LEPTON MIXING: 

We argued above that lepton mixing is a consequence of an intermediate meson/neutrino 
state which couples to both charged leptons. The intermediate meson state is a low energy 
or "long distance" representation of the coupling of charged leptons to quarks via charged 
current interactions and is akin to the mixing between K K via intermediate states with 
two and three pions. This "long distance" (low energy) contribution to the charged lepton 
self energy is depicted in fig. (JTJ). 

M~ 

'•••••>.••••*' 

v 3 

FIG. 1: Long distance contribution: intermediate state with Uj and M = ir, K. 

This is a low energy representation of physical process in which a lepton couples to an 
intermediate W vector boson and a neutrino mass eigenstate, followed by the decay of the 
(off-shell) W into quark-antiquark pairs with the quantum numbers of the pseudoscalar 
mesons. Therefore we also expect a short distance contribution in which the intermediate 
state corresponds simply to the exchange of a W boson and a neutrino mass eigenstate. This 
contribution to the charged lepton self-energy is depicted in fig. (J2]). 




FIG. 2: Short distance contribution: intermediate state with Vj and W . 

We shall calculate both self-energy diagrams to properly ascertain each contribution to 
the full charged lepton propagator. 

However, before carrying out the detailed calculations, we point out that there are also 
electromagnetic and neutral current contributions to the self-energies. However these are 
flavor diagonal, thus while they will both contribute to the self-energies, only the charged 
current contributions ( long and short distance) lead to off diagonal self energies which lead 
to charged lepton mixing. Furthermore, both long and short distance self-energies are of the 
general form 

Vafi OC Y^UajSjUjp (IH.l) 
3 
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where Sj is the propagator of neutrino mass eigenstates. Unitarity of the mixing matrix, 
J2j UajUjp = 5 a p, leads to a GIM (Glashow-Ilioupoulos-Maiani) type-suppression of the off- 
diagonal matrix elements: if all the neutrinos in the loop are degenerate, unitarity entails 
that there is no off-diagonal contribution to the self-energy; furthermore, this argument also 
suggests that the off-diagonal terms will be dominated by the most massive neutrino state. 



A. Short distance contribution: 

We begin by computing the self energy contribution from W exchange depicted in fig. 
d2J). Throughout this calculation, we shall be working in the physical unitary gauge and in 
dimensional regularization. Upon passing to the basis of mass eigenstates that define the 
neutrino propagators, ip Ua = J2jU a jipj, the charged current contribution to the self energy 
matrix is given by 



a/3 



V2 



E 



d 4 k 



£(# + rrij) 
k 2 — m 2 + it 



, . n fiv _ 
H9 



(p-k) 



2 - M 2 
lvl w 



(III.2) 



where L, R = (1 =p 7 5 )/2 respectively. 

This integral is calculated in dimensional regularization. We introduce a renormalization 
scale k, which we choose k = My/ thus renormalizing at the W — pole, and define 



= -p 2 x(l -x) + m 2 x + M^(l - x) , 
separating explicitly the divergent and finite parts in the MS scheme we find 
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I?(p 2 ;x)+I (p 2 ;x) 



dx 



(III.3) 



(III.4) 



where 



I-(p 2 ;x) 



g 2 M, 



w 



2(47r) J 



X) 



3A 7 - , 

2+ M^ + (1 
lvl w 



X 



p 



M 2 W1 



7 + ln47r 



g 2 M } 



w 



2(4tt) 



x 



2 + 



3Aj 



+ d 



X 



M 2 W 



Aj 2xm 2 
In -r4r + 



M 2 W 



M 2 W1 



(III.5) 



(III.6) 



Unitarity of the neutrino mixing matrix in the form ^ . U a jU*p = 5 a p leads to GIM-like 
cancellations in both the divergent and the finite parts for the off-diagonal components of 
the self-energy matrix. Therefore, for a ^ (3, we find 



S a/8(P) 



where 



g 2 M, 



w 



64vr 2 



5>a;^ 



IV1 W 



7 + ln47r) ; a ^ f3 



(III.7) 



(III.8) 
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and 



9 2 M } 



w 



327T 2 




2 + —^ + 



2m) 



x > (1 — x) dx ; a ^ (5 



xf 



p 



M 2 W 



In 



M 2 W 



+ 



(III.9) 



B. Long distance contribution: 



We now turn our attention to the intermediate state described by the exchange of a 
7i /K meson and a neutrino mass eigenstate. This is the state that suggested charged lepton 
mixing from the density matrix treatment in the previous section. A difficulty arises in the 
calculation of the meson exchange because in order to properly describe the coupling between 
the meson and the charged lepton and neutrinos, we would need the full off-shell form factor 
Fm{q 2 ) which is a function of the loop momentum since the meson is propagating off its 
mass shell in the intermediate state. Clearly this is very difficult to include in a reliable 
calculation, therefore we restrict our study to an estimate of this contribution obtained by 
simply using the on-shell value of the form factor, namely the meson decay constant F M in 
order to obtain an admittedly rough assessment of its order of magnitude. 

Under this approximation the contribution to the self-energy matrix from this interme- 
diate state is given by 



d A k 



+ rrij) 
k 2 — m 2 + ie 



(p - k) 2 -Mlj + ieJ 
(111.10) 

where Mm is the meson mass. The width of the meson may be incorporated via a Breit- 
Wigner approximation M M — > M M — iT M /2, however this will only yield a contribution 
which is higher order in Gf- 

The calculation is performed in dimensional regularization, choosing the renormalization 
scale k = Mw as for the short distance contribution, introducing 



S j = -p 2 x (l -x) + m)x + M M (1 - x) , 



(111.11) 



and separating the divergent and finite parts in the MS scheme we find 



Y. U *i U h / Uf(p 2 ;x) + Jj(p 2 ;x) 



dx 



(111.12) 



where 



Jf(p 2 ;x) 



M 2 w e F 2 M 
(4tt) 2 



J-(p 2 ;x) 



M 2 w e F 2 M 



(4vr) s 



2x 



2 u 3 



^■(l + 3x) - (1 -x)x 2 p 
6< 



M 2 W 



+ 



x)x 



P 
M 2 W 



7 + In 47r 



;i + 3x) Sj 



M 2 W 



In 



M 2 W 



(111.13) 
(111.14) 
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For the off-diagonal matrix elements, unitarity of the neutrino mixing matrix leads to GIM 
type cancellations as in the short distance case, therefore for a 7^ (3 we find 



^(p) = ^[i + i(p 2 )] 



where 



6 32n 2 



1VI W 



7 + In An 

e 



M 2 w ~ e F 2 M 

167T 2 







2x 



! U J I 1 1 \ 2 P 

M? 



1 w 



H (l-x)x 2 



;i+3x) 



a ^ /9 



In 



M 2 W 



However, with Fm oc GV/t^k" and ~ 100 MeV it follows that 

2 



/ -2 1 r 2 2 / 9fn,k 

F M M W <xg 



Mi 



~ 10~V 



(111.15) 

(111.16) 

<ix ; a ^ (3 
(111.17) 

(111.18) 



therefore the long distance contribution is negligible as compared to the short distance 
contribution and, to leading order, the off-diagonal components of the self energy are given 
byeqns. ([ITO- MM) ■ 

As noted previously, unitarity of the neutrino mixing matrix entails that the flavor off 
diagonal matrix elements of the self-energy vanish either for vanishing or degenerate neutrino 
masses. Obviously the contribution from light active-like neutrinos is strongly suppressed by 
the ratios m 2 /M^, hence these off-diagonal matrix elements are dominated by the heaviest 
species of sterile neutrinos. 

Thus charged lepton mixing is enhanced by intermediate states with heavy sterile neutri- 
nos. This is one of the main results of this article. 

// even the heaviest generation of sterile neutrinos feature masses mj <C My/ and for 
p 2 <C My? the following order of magnitude for the off-diagonal component z^ e is obtained 



2 

0"°j ' 



(HI.19) 



as it will be seen below this estimate determines the mixing angles up to kinematic factors. 



IV. FULL PROPAGATOR: MIXING ANGLES AND PROPAGATING MODES. 
A. Full propagator and mixing angles: 

To treat //, e mixing it is convenient to introduce a flavor doublet 
The general structure of the self-energy is of the form 



S(p) = z L (p 2 )p + 5M L (p 2 



L + 



i R (p 2 )f + SM R (p 2 



R. 



(IV.l) 



(IV.2) 
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The neutral current interactions contributes generally to the right and left components of 
the self-energy but are diagonal in flavor and as are the electromagnetic contributions. The 
V — A nature of the charged current interactions is such that their contribution is only of 
the form z^(]9 2 )^L and is the only contribution that yields flavor off-diagonal terms and are 
ultimately responsible for // — e mixing. To cancel the poles in e in the self-energy we allow 
counterterms in the bare Lagrangian 



C ct = ^(5Z ct - l)pV + ^5M^ + h.c. . 
The full propagator S now becomes a 2 x 2 matrix which is the solution of 

|^l+^Z ct -l)-E(p)-M]S = l 
where the boldfaced quantities are 2x2 matrices and 



M 



M M 
M e 



(IV.3) 



(IV.4) 



(IV.5) 



In what follows we will assume that that M contains the renormalized masses and we will 
neglect finite momentum dependent contributions to M since these will only generate higher 
order corrections to the mixing matrix, as will become clear below. 

We will choose the counterterm (SZ ct — 1) in the MS scheme to cancel the term z^r i n 
eqn. (JTTTT7J) . Therefore equation (1IV.4|) becomes 



where 



■Z^R-M 



J L,R 



(IV.6) 



(IV.7) 



The leading contribution to the off-diagonal matrix elements is given by the "short-distance" 
term eqn. (1IH.9|) . 

Multiplying on the left both sides of flTV6]) by f + MZ fi L + MZ L R and writing the full 
propagator as 

S = RS R + L S l 



where 



we find 



s R = A R (p 2 )\p + B R (p 2 : 

S L = A L (p 2 )^ + B L (p 2 ) 
[p 2 Zr 1 ~ MZ L Mj A/jQo 2 ) 

p 2 Zl 1 -mz r m)a l (p 2 ) 



(IV.8) 

(IV.9) 
(IV. 10) 

(IV.ll) 
(IV.12) 



and the conditions 



B R (p 2 ) = MZ L (p 2 ) ; B L (p 2 )=MZ R (p 2 



(IV. 13) 
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In what follows we will neglect CP violating phases in U a j with the purpose of studying 
fx — e mixing in the simplest case. Under these approximations we find 

I): The solution for A R (p 2 ) in eqn. ( llV.llj) is obtained as follows. Consider the diago- 
nalization of the inverse propagator 



p 2 Z^-MZiM 



where 



1 



Q r (p 2 )+Q r (p 2 ) 



X R {p 2 ) fcos2 9 R {p 2 ) sin2 9 R {p 2 ) 
sin2 9 R (p 2 ) - cos 2 9 R (p 2 



Q> 2 " 
QfiP 



mi 



P 2 -M 2 U [Z R 



e 



Zn 1 



p 2 -M 2 [Z R ) (Z, 



ee \ / ee 



and 



{Q*{p 2 ) - Q*{P 2 ))' +±(M,M e z[^p 2 



To leading order we find the mixing angle to be given by 

tan 2 9 R {p 



., _ 2M,M e z{^ e (p 2 ) 



,/ze * 

M 2 - M* 



The matrix above can be diagonalized by a unitary transformation 

U[9} = 

in terms of the mixing angle R (p 2 ), namely 



cos 9 sin 9 
- sin 9 cos 9 



U[9 R {p 2 )} 
where 



p 2 Zp 1 - MZ L M 



U- l [9 R {p 2 )] 



Q R (p 2 ) - qr{p 2 ) 









QfiP 2 ) + MP 2 



g R (p 2 ) = (M 2 -M 2 )sm 2 29 R (p 2 ) 



and to leading order we find 



A R {p 2 )=U- 1 [9 R {p 2 )] I 9W-teG*H* 
which (to leading order) simplifies to 

A R (p 2 ) -U-^Rip 2 )} P J -^(P 2 )-^)+ie 





1 



U[9 R {p 2 )} 



Q§{p 2 )+en{p 2 )+it 







p 2 -Mi(p 2 )+g R (p 2 )+ie 



(IV.14) 

(IV.15) 
(IV.16) 

(IV. 17) 



(IV. 18) 



(IV.19) 



(IV.20) 
(IV.21) 

(IV.22) 



U[9 R {p 2 )} . (IV.23) 



In the above expressions M 2 (p 2 ), M 2 (p 2 ) include the finite renormalization from the diagonal 
contributions of the self-energy matrix which have not been calculated here, furthermore 
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the residues at the poles (wave-function renormalization) are also finite since the (local) 
divergent contributions are canceled by the counterterm. 

Therefore Sr can be written in the basis that diagonalizes the kinetic term 



p 2 -M 2 {p 2 )-Q R {p 2 )+ie p 2 _ M 2 ( p 2 ) _ — ( p 2 ) + ie 
\ p 2 -M 2 (p 2 )+ SR (p 2 )+ie p 2 -M 2 (p 2 )+g R (p 2 )+it / 



where 



b R (p 2 ) = U{6 R (p 2 )] MZ L (p 2 )U- l [9 R (p 2 )} 



(IV.24) 



(IV.25) 



II): We proceed in the same manner for Al(p 2 ), namely consider diagonalizing the inverse 
propagator 



p 2 Z- 1 - MZ,M = - [Q> 2 ) + Qe(p ) 
where 



Xl(p 2 ) ( cos2 9 L (p 2 ) sin2^(p 2 ) 
sin2# L (p 2 ) -cos2# L (p 2 



Qe(p 2 ) : 



/'/' 



Zl 1 



p 2 -M 2 [Z H 
p 2 -M 2 [Z u 



/'/' 



and 



Hp 2 ) 



(Qi(p 2 )-Q L e(p 2 )) 2 + ^P 2 4,,e(p 2 



Again, to leading order we find the mixing angle to be given by 

_ 2p 2 z{ ue {p 2 ) 
tan 2 9 L (p 2 ) - ^ 1 



M?-M? ' 



(IV.26) 

(IV.27) 
(IV.28) 

(IV.29) 
(IV.30) 



The matrix above can be diagonalized by the unitary transformation ( 1IV.19j) now in terms 
of the mixing angle 6l{p 2 ), namely 



U[6 L (p 2 )] [p 2 Z L X -MZrMJw- 1 ^ 2 )] 
where 



Q l ,(p 2 ) - ql{p 2 







QeiP 2 ) + QL{P 2 



Q L (p 2 ) = (M 2 -M 2 ) sin 2 29 L (p 2 ) 



and to leading order we find 
A L {p 2 )=U- l [6 L {p 2 )} 



i 



}jl(P 2 )-eL(p 2 )+ie 






1 



U[9 L {p 2 )], 



(IV.31) 
(IV.32) 

(IV.33) 



Qi{P 2 )+8L(p 2 )+ie 



Neglecting the diagonal contributions to mass renormalization, but keeping the (finite) 
wave function renormalizations, the result (lIV.34j) simplifies to 



zUp 2 ) 



Al^^U-^O^p 2 )} P 2 -^iP 2 )sL(P 2 )+ie 





zL(p 2 ) 

p 2 -M 2 {p 2 )+g L (p 2 )+ie 



U[9 L (p 2 )), (IV.34) 
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Just as in the previous case, M 2 (p 2 ), M 2 (p 2 ) include the finite contribution from mass 
terms in the self energy and the residues at the poles are also finite, the local, divergent 
contribution being canceled by the counterterm. 

The component can now be written as 



/ 



U[9 L (p 2 )]S L ir L [6 L (p')) 



^> 2 )^ e (P 2 ) 



\ 



p 2 -M 2 (p 2 )-g L (p 2 )+ie p 2 - M 2 (p 2 ) - g L (p 2 ) + ie 

ZL(P 2 ) [p+bL(p 2 )] 

p 2 -M 2 {p 2 )+Q L {p 2 )+ie / 



zUp 2 )*%kp 2 ) 



\ p 2 -M 2 (p 2 )+g L (p 2 )+i 



(IV.35) 



where 

b L (p 2 ) = U[9 L {p 2 )} MZ R {p 2 ) U- x [9 L {p 2 )\ . (IV.36) 

An important aspect is that the mixing angles n{p 2 ) , @ l{p 2 ) not only are different for 
the R, L components a consequence of the V — A nature of charged currents, but also that 
they feature very different momentum dependence, 



On(p 2 ) - ^H> 2 ) 



Ml 



(IV.37) 



Near the muon mass shell p 2 ~ it follows that 9l ^> Or, for near the electron mass 
shell p 2 ~ M 2 it follows that 9 R 3> 9 L . Off-shell, for virtuality p 2 3> mixing of the L 
component becomes dominant. 

In general, the transformations that diagonalize the kinetic terms f for both the positive 
and negative chirality components do not diagonalize the mass terms. In the basis in which 
the kinetic terms are diagonal, the pole-structure of the propagator is revealed and the 
propagating modes can be read-off. This basis, however, does not diagonalize the mass 
term of the propagator and attempting to diagonalize the latter either via a unitary or a 
bi-unitary transformation will lead to an off diagonal matrix multiplying the kinetic term. 
A similar situation has been found in different contexts 



B. Propagating modes: the effective Dirac equation: 

The nature of the propagating modes is best illuminated by solving the effective Dirac 
equation for the flavor doublet, which corresponds to the zeroes of the inverse propagator, 
namely 

f Zl 1 L + f Z B X R - M (p) = . (IV.38) 

It is convenient to work in the chiral representation and expand the positive and negative 
chirality components in the helicity basis 

B ■ ^rv h (p) = hv h {p) ; h = ±1 (IV.39) 

Ipi 

in terms of which the spinor flavor doublet \1/ becomes 

*(p) = XS® (If ) ' ( IV - 4 °) 
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where £ h ' are flavor doublets that obey the following equations 

Zl\p + ph)ti + Mtf = (IV.41) 

Z^(po + V h)tf + = . (IV.42) 

The positive and negative energy and helicity components are given by (p = \p\) 

( z£m^J ; p >0;/i = l ; ( ^ ) ' Po>0,/i=-l (IV.43) 

J ; R,<0;/i = l ; (^zjfl^j ; p >0,h = -l. (IV.44) 
The flavor doublets obey 

(/Z^-MZiM)^^) = (IV.45) 

(p 2 7,-j} - MZ fl M)( L (p) = , (IV.46) 

using (IIV.201 IIV.31[) we find that the rotated doublets 



obey the following equation 

V o q? (p 2 ) + te (p 2 ) ) \<p*(p) ) -° (IV - 48) 



Q l ,{p 2 ) - ql{p 2 ) o W*) 



. (IV.49) 



Neglecting perturbative renormalization of the /i, e masses, for p 2 ~ M 2 the propagating 

modes correspond to <p x :,L ^ ; ^ L = and the mixing angles for R, L components are 
6 RL (M 2 ) respectively, with 

M*^4W ; Ol{K)-£AK) ( iv - 5 °) 

defining the /i-like propagating modes 

&) \ - ^l (v) ( cose L (M 2 ) \ . / C?(P) A _ *, x ( cos^(M 2 ) \ 

Similarly for p 2 ~ M 2 the propagating modes near the electron mass shell correspond to 
ip 2 ' 7^ ; Lp x ' =0 and the mixing angles for R, L components are 9r,l{MV) respectively, 
with 

M M 2 
Q R {Ml)^^z{^Ml) ; 9 L {Ml) ~ j^zl^Ml) (IV.52) 
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defining the relation between the flavor doublets and the propagating modes on the respective 
mass shells, namely 



#(P) 



-sm8 L {M 2 
cos0 L (M 2 ) 



-sin^(M 2 
cos^(M 2 ) 



(IV.53) 



The expressions (1IV.51IIV.53P combined with (lIV.43|IV.44p give a complete description of 
the propagating modes. 
Discussion: 

There are several noteworthy aspects of the above results: 

• The unitary transformation that diagonalizes the kinetic term i> in the propagator does 
not diagonalize the mass term. This can also be seen in the spinor solutions because 
the terms M in the spinors (IIV.43piV.44j) are generally not diagonalized by the 
rotation matrices U[9 LjR ]. 

Furthermore, 



A similar situation has been found in different contexts in refs. 73, 74 



the momentum dependence of the mixing angles entails that the relation between 
flavor doublets and the propagating modes is not a simple unitary transformation but 
requires different mixing angles for different chiralities and mass shells as highlighted 
byeqns. (IIV.5imV.53p . 



The mixing angles are GIM suppressed favoring heavier neutrinos in the intermediate 
state and are momentum and chirality dependent. This means that off-shell processes 
necessarily mix charged leptons with virtuality and chirality dependent mixing angles. 
For p 2 -C Myy and assuming that the heaviest sterile neutrinos feature masses <C M^, 
from eqn. ( 1III.19P we find the positive and negative chirality mixing angles for \i — e 
mixing 



G F M e \ - 



4tt 2 M„ 



0l( P 2 



4vr 2 M 2 

(J, 



E w 



(IV.54) 



thus the mixing angles are dominated by the heaviest generation of neutrinos. In 
particular if heavy sterile neutrinos do exist, these new degrees of freedom will yield 
the largest contribution to charged lepton mixing. 

Considering, for example, that there is only one generation of heavy sterile neutrinos 
with mass Ms, and assuming that ~ U e i, the recent results from the PIENU 
collaboration at TRIUMF 56J reporting an upper limit \U e i\ 2 < 10 -8 (90% C.L.) in 



the neutrino mass region 60 — 129MeV/c 2 allows us to estimate an upper bound for 
the negative chirality mixing angle near the /x mass shell, 



6 AM 2 ) < 10" 14 — 

V ^' ~ VlOOMeV/ 



Mc 



(IV.55) 



V. RELATION TO LEPTON FLAVOR VIOLATING PROCESSES: 

Charged lepton mixing via intermediate states of charged vector bosons and neutrino 
mass eigenstates are directly related to lepton flavor violating processes. An important 
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process that is currently the focus of experimental searches 67, 68] and a recent proposal [6 



is the decay // — > which is mediated by neutrino mass eigenstates [75l-l78l| ; the importance 



of heavy sterile neutrinos in this process has been highlighted in ref . [79( . However, to the 
best of our knowledge, the relationship between this process and a mixed \i — e propagator 
has not yet been explored. Such relationship is best understood in terms of the three-loop 
muon self-energy diagram in fig. fl2Ka)), the Cutkosky cut along the intermediate state of 
the electron and photon yields the imaginary part of the muon propagator on its mass shell, 
and determines the decay rate fi — > cy, this is depicted in fig. (j3]-(b)). 




FIG. 3: Lepton flavor violation: fig. (a): three loop contribution to the Cutkosky cut through 
the photon and electron intermediate state yields the imaginary part describing the flavor violating 
decay u — > ej of fig.(b). 



However, the self-energy diagram ([3]-(a)) is only one diagonal component of the full // — e 
self-energy, the corresponding three loop diagram for the off-diagonal component is shown 
in fig. H 




FIG. 4: Three loop contribution to which is the off-diagonal counterpart of fig. ([3^a). 



Because of the different external particles, a Cutkosky cut of this diagram through the 
photon and electron internal lines cannot be interpreted as a decay rate. However, this anal- 
ysis clearly indicates the relationship between /i — > cy, a distinct indicator of lepton flavor 
violation, and charged lepton mixing in self-energy diagrams, both a direct consequence of 
neutrino mixing. 

We note that whereas the branching ratio for \i — \ ej is oc GVa| J2j U[_ijUj e \ 2 m?j we find 
that the one- loop mixing angles are momentum dependent, different for different chiralities 
and the largest angle for on-shell states corresponds to the negative chirality muon-like 
combination, in which case the angle is of order Gf Ylj U^U^m 2 . 
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Possible other contributions: The diagram in fig. ([3]-(b)) suggests that \i — e mixing 
may lead to further contributions. Consider the flavor blind electromagnetic vertices of \i 
and e, if the mixing angles were momentum independent, unitarity of the transformation 
would entail a GIM cancellation between off-diagonal terms in the electromagnetic vertices, 
just as for neutral currents. However, muon-like and electron-like mass shells feature very 
different mixing angles which suggests that off diagonal contributions arising from replacing 
the /i and e fields in the electromagnetic vertices by the correct propagating states would not 
cancel out because of different mixing angles. This can be seen from the relation between 
the propagating states and the /x, e states given by eqns. flIV.51IIIV.53l) . writing 



ipn = cos #1^1 - sin8 2 ip 2 ; ip e = cos 6 2 ip 2 + sin^i^i (V.l) 

respectively for positive and negative chirality components with the respective angles On = 
9l(M 2 ); 9 2 l = Ol(M 2 ) etc., it follows that the electromagnetic vertices feature a mixed term 
of the form 

oc Tp 2L ^A^ 1L {B lL - 9 2L ) + L^R. (V.2) 

If the mixing angle (s) were momentum independent Q\ = 9 2 and this term would vanish in 
a manner similar to the GIM mechanism. Therefore mixing with momentum and chirality 
dependent mixing angles suggests that the contribution to /i — >■ e~f from the vertex (IV. 21) 
depicted in fig. (JSJ) becomes possible. 



FIG. 5: Further contribution to S^e from \i — e mixing. 

This contribution differs from that of fig. (j3]-(b)) in two major aspects: i) rather than 
an extra W propagator, it features an electron propagator in the intermediate state, which 
would suggest a large enhancement with respect to the usual contribution, ii) a very small 
mixing angle which suppresses the enhancement from the electron propagator in the inter- 
mediate state. Thus a detail study of both effects and their impact is required for a firmer 
assessment. 

This argument, however, needs to be scrutinized further by analyzing the imaginary part 
of the propagators keeping both the diagonal electromagnetic contribution as well as the off 
diagonal charged current contribution. Upon the diagonalization of the propagator there is 
an interference between the diagonal and the off-diagonal terms, this can be gleaned from 
(lIV.14pV.26p . The imaginary part of the propagator evaluated at the mass shell of the 
muon-like propagating mode, namely p 2 ~ M 2 , would yield the contribution to the diagram 
of fig. (E]) from the interference between the diagonal electromagnetic contribution which 
features an imaginary part for p 2 > M 2 and the off-diagonal charged current contribution. 
Since this is a contribution to the self energy of higher order than the ones considered 
here, a firmer assessment of this new contribution merits further study and will be reported 
elsewhere. 
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VI. CONCLUSIONS AND FURTHER QUESTIONS 



In this article we studied charged lepton oscillations and mixing. The decay of pseu- 
doscalar mesons leads to an entangled quantum state of neutrinos and charged leptons (we 
focused on 7r, K decay leading to /i, e). If the neutrinos are not observed, tracing over their 
degrees of freedom leads to a density matrix for the charged leptons whose off-diagonal 
elements in the flavor basis reveals charged lepton oscillations. While these oscillations de- 
cohere on unobservably small time scales < 10 -23 s, we recognize that they originate in a 
common set of intermediate states for the charged leptons. This realization motivated us 
to study the mixed \i — e self-energies and we recognized that charged- current interactions 
lead to a dominant "short distance" contribution to \i — e mixing via W-exchange and an in- 
termediate neutrino mass eigenstate, and a subdominant (by a large factor) "long distance" 
contribution to mixing via an intermediate state with a pseudoscalar meson and neutrino 
mass eigenstate. We include the leading contribution in the propagator matrix for the \i — e 
system focusing on the off-diagonal terms which imply \i — e mixing. We find that the mixing 
angles are chirality and momentum dependent, the chirality dependence is a consequence of 
V — A charge current interactions. The kinetic and mass terms of the mixed propagator 
cannot be diagonalized simultaneously either by unitary or bi-unitary transformations as 
they feature different matrices and spinorial structure. Diagonalizing the kinetic term dis- 
plays the poles which describe muon-like and electron-like propagating modes for which we 
find explicitly the wave functions, but the mixing angles evaluated on the respective mass 
shells (and chiralities) are very different. We find the positive and negative chirality mixing 
angles for p 2 <C M^ to be approximately given by 

M 3 M 3 

therefore dominated by the heaviest generation of sterile neutrinos. For one (dominant) 
generation of massive sterile neutrinos with mass Ms, the recent results from the PIENU 
collaboration at TRIUMF{56[, suggests 

We discussed the relationship between the lepton flavor violating decay // —> ej, the focus 
of current searches (oTI. 68| and proposalsj69|, and charged lepton mixing, pointing out that 



a positive measurement of the former confirms the latter. Furthermore, we advance the 
possibility of further contributions to \i — » ey arising from the fact that the \l — e mixing 
angle is momentum dependent and differs substantially on the mass shells of the propagating 
modes voiding a GIM mechanism for the electromagnetic vertices. 

The (four) momentum dependence of the /x — e mixing angle may be the source of novel 
off-shell effects whose potential observational manifestation merits further study. We expect 
to report on ongoing study on these issues elsewhere. 

Furthermore, in order to present the main arguments in the simplest case, in this article 
we have not considered CP- violating phases in the mixing matrix elements U a j, including 
these phases merit further study since this aspect could indicate potentially rich CP-violating 
phenomena from the charged lepton sector induced by CP- violation from the neutrino sector 
which merits further and deeper study. 
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